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• Abstract This expository review is devoted to fish swimming and bird/insect flight, (i) The simple 
' waving motion of an elongated flexible ribbon plate of constant width propagating a wave distally 
. down the plate to swim forward in a fluid, initially at rest, is flrst considered to provide a fundamental 
concept on energy conservation. It is generalized to include variations in body width and thickness, 
with appended dorsal, ventral and caudal flns shedding vortices to closely simulate flsh swimming, for 
Q> ■ which a nonlinear theory is presented for large-amplitude propulsion, (ii) For bird flight, the pioneering 
studies on oscillating rigid wings are discussed with delineating a fully nonlinear unsteady theory for a 
2D flexible wing with arbitrary variations in shape and trajectory to provide a comparative study with 
experiments, (iii) For insect flight, recent advances are reviewed by items on aerodynamic theory and 
modeling, computational methods, and experiments, for forward and hovering flights with producing 
leading-edge vortex to yield unsteady high lift, (iv) Prospects are explored on extracting prevailing 
intrinsic flow energy by fish and bird to enhance thrust for propulsion, (v) The mechanical and biological 
principles are drawn together for unified studies on the energetics in deriving metabolic power for animal 
locomotion, leading to the surprising discovery that the hydrodynamic viscous drag acting on swimming 
■ fish is largely associated with laminar boundary layers, thus drawing valid and sound evidences for a 
resounding resolution to the long-standing fish-swim paradox proclaimed by Gray (1936, 1968). 

Sl:^ ■ 1. INTRODUCTION 

Fish swimming and bird flight are of agelong human interests, with their fascinating nature to 
' appreciate, diversity to behold, and complexities to comprehend. The focal inspiration for the interest 
, is curiosity. That is the seed of creativity, thereby bringing forth fruitful formation of sound physical 
conceptions, successful developments of required mathematical tools, and novel engineering. Without 
curiosity, one may wonder what would have become to riding in jet airplanes and in cruising ships across 
oceans for the ease and comfort in transportation. This is how the advances in science and technology 
r> ' have been all about. 

5^ I Oscillatory airfoil theory has been developed since 1920s by such pioneers as Wagner (1925), 

Theodorsen (1935), Kiissner (1935) and others, given in one or another form of linear theory. The 
application was flrst to aero-elasticity problems to expound the fluttering forces and their effects on 
airplane wing-fuselage system flying through turbulent and gusty air. The fundamental elements of the 
theory would also be applicable to studies on flsh swimming and bird flight. 

Keen interests have been stimulated by the paper published by G.I. Taylor (1952) on the swimming 
of long and narrow animals, arousing research activities at various institutions, especially at Cambridge 
University and Caltech. Close academic interactions and collaborations between various centers plus 
their prolific publications rendered it timely to hold the first "International Symposium on Swimming 
and Flying in Nature" in Aug. 1974 at Caltech, participated by 250 of the active workers in the field. In 
Sept. 1975, another symposium, called by Lighthill and Weis-Fogh, was held at Cambridge University 
on "Scale Effects in Animal Locomotion" , a new subject of great significance. These symposia are of 
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value to have set fruitful directions for vital advances in this growing important field. 

For fish propulsion, the first conspicuous feature is the undulatory motion of the flexible fish body 
passing a wave distally in swimming. This feature is abstracted to expound the simplified geometry of 
an elongated flexible ribbon plate of constant width, to which the slender-body analysis is applicable as 
presented in §2.1-2.3. The fundamental quantities derived stepwise in this simple case provide a crisp 
clear physical concept on momentum, power expended for generating thrust on the plate with some 
energy expelled at an optimum efficiency for the plate to swim forward. These primary results are then 
generalized in §2.4.1-2.4.3, by including all the other main factors to simulate fish swimming in nature 
on general slender-body theory given in §2.5. To this end, the nonlinear theory introduced by Lighthill 
(1971) for large-amplitude swimming is finally developed in this work to completion in §2.6. 

For bird flight, the classical oscillating airfoil theory is discussed in §3.1, followed by presenting a 
fully nonlinear unsteady theory developed by Wn (2001b)- (2007) for a flexible two-dimensional wing 
with arbitrary variations in wing shape and its trajectory as summarized in §3.2, thus furnishing a 
comparative study by comparisons with two well noted experiments on heaving and pitching of airfoils, 
respectively. The insect flight is reviewed in §4, itemized for the aerodynamic theory and modeling in 
§4.1, computational fluid dynamics in §4.2, and experimental investigations in §4.3. 

For fish swimming and bird flight in wild nature, e.g. in torrential rivers and in whirling winds, 
the prospects are explored in §5 on the animal's apt capabilities of extracting intrinsic flow energy to 
gain more thrust for propulsion in unfavorable environments with aspects relevant to control theory. 

For the energetics of animal locomotion examined in §6, the most effective and efficient parts of 
mechanics (on reactive theory for determining thrust) and biology (on measuring the metabolic rate 
of converting oxygen consumption into muscular power) are drawn together for unified studies on the 
total energy conservation and the balance between the thrust exerted by the fish and the hydrodynamic 
viscous drag acting on the fish. The highly rewarding results are as just summed in Abstract. 

2. FISH SWIMMING 

In general, a fish swimming in water adopts a sidewise undulatory waving motion of its flexible 
body of length i and depth 2b say, commonly appearing in the form of the so-called slender body, or 
elongated body, typifled by the slenderness-ratio S = b/i <^ 1. The waving motion in general propagates 
a wave distally down the body to generate a forward thrust T to balance out the hydrodynamic viscous 
drag D acting on the fish so as to sustain a steady swimming speed U as so resulted. For most differing 
species of grown-up fishes in cruising swim, the Reynolds number Re = Ul/u {v being the kinematic 
viscosity of the fluid) would be so high (~ 10"^ for small slow fish, ~ 10^ for whales) and the boundary 
layer enveloping the fish so thin that its thickness can be neglected for evaluation of thrust T and power 
P of fish locomotion, a powerful method known as the reactive theory of animal locomotion, which is 
commonly adopted with success. With body motion specified and analyzed with certain viscous drag 
distributed along the body, it may lead to found a resistive theory of animal locomotion, of which the 
basic ideas have been delineated by Lighthill (1973a, §2). On the other hand, the bio-chemical studies of 
the oxygen consumption by a steadily swimming fish and its metabolic conversion into muscular energy 
for swimming is a science and technology well achieved, as shown in their precious contributions to this 
problem. Perhaps due to this contrast in recognizing their early status quo, a highly inspiring paradox 
on fish locomotion was proclaimed by Sir James Gray (1936, 1968). Briefly, it asserts that either fish 
can experience a viscous drag much smaller than the drag acting on their artificial smooth model or they 
can deliver muscular power in sustained swimming an order of magnitude greater than that delivered by 
the same mass of muscle of warm-blooded animals working at the same level of activity. This claim has 
provided a strong stimulus to have attracted both mechanical and biological studies on fish locomotion. 
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2.1. Simplified Reactive Swimming 

In aquatic animal locomotion, the first conspicuous feature to be noted is the flexibility of fish body. 
An early study on this feature is a Fourier synthesis of a waving plate pursued by Wu (1961), with plate 
undulation represented by a Fourier series in general. In addition, another feature of acceleration of 
waving plate was expounded by Wu (1962). As explained, these general geometric shape functions 
are first abstracted here to a generic undulatory plate of a very simple geometry for fundamental 
study. Thus, we consider the propulsion, on reactive theory, of a flexible ribbon plate of constant width 
spanning ~b < y < b, initially stretched straight from the leading edge at x = to the trailing edge 
aX X = £ {6 = b/i <^ 1), with its thickness neglected. It then starts from time t = to propagate a 
continuous waving motion given in cartesian coordinates x = (x, y, z) flxed in the body frame as 

z = h{x,t) {0<x<i, \y\<b, t>0), (1) 

with the plate immersed in a fluid of density p, which is unbounded and has a uniform free stream 
velocity {U, 0, 0) far away. Here, U is the swimming velocity of the plate in the — x-direction in the 
absolute frame with fluid at rest far away. The fluid is assumed incompressible and inviscid, and the 
resulting flow assumed irrotational except for the vortex sheet generated at the ribbon surface and the 
vortex sheet shed from the trailing edge to become free vortices drifting away with the local fluid. 

For flow analysis, we adopt the 'linear slender-body theory^ applicable to bodies of S = b/£ <^ 1 
moving with gradual gradients so that the fluid motion at arbitrary station x is all determined in the 
yz-plane transverse to the x-direction. Thus, the flow at a given x at time t depends only on its own 
h{x,t) and otherwise unaffected by the plate motion both up- and down-stream. 

Thus, originally, the 3D flow has a perturbation velocity potential (p{x, y, z) such that 

= {u,v,w), q = {U + u,v,w), = 0, and ^ — )■ as |a;| — > oo, (2a) 

where q is the total velocity vector, and the Bernoulli equation for pressure p is given by 

p/p = _Dcl>-^{V<t>f, D = ^^ + U^. (2b) 

However, for the linearized slender-body flow in the yz-plane, or in the complex C = (y + «-2)-plane, 
there exist a complex velocity potential f = (p + iij: and complex velocity x = df /dC = v — iw, or 

C = y + iz, f = (j) + iil), x = df/dC = V - iw; Vy + Wz = 0, V;,-Wy = 0, (3) 

where f and x are analytical functions of one another under the incompressibility and irrotationality 
conditions as expressed by the final two relations in (3), where Vy = dyV = dv/dy, etc. 

The plate motion (1) induces a flow velocity field with its z-component at the plate given by 

w±=w{x,y,h±0,t) = Dh{x,t) = {dt + Ud^)h = W{x,t) (!y| < 6). (4) 
With this boundary condition at the plate and with % — ^ as |C| — ^ oo, the fiow has the solution: 

^ = y-inj = iw(^-^^==-l^ , f = (P + ii^ = iWiVe -b'^-O {cP± = TWVb'^-y^). (5) 

The flow momentum has zero y-component since v is odd in y, and its z-component per unit longitudinal 
length (or unit thickness in x) is given by 

M = pJ J ^dydz = pj{(j)-- (p+)dy = 2pW J ^b"^ - y^dy = mW{x, t) (m = pirb"^), (6) 
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where the 5-integral is over the entire yz-plane per unit length in x, giving, by the z-integral, the jump 
{(p- — (/)-!-) which exists only over the plate, hence the result, with m = pnb'^ signifying the added mass 
(or the virtual mass) of the fluid at station x. In terms of flow singularities, this momentum is due to 
an elliptic distribution of doublets (or dipoles) at the plate, namely {(j)- ~ 4'+) = ^Wy^b"^ — y^. 

The instantaneous lift per unit length of the plate is the force, positive in the 2:-direction, called 
the specific lift, C{x,t), given by 

/b I'b 
{p- - p+)dy = -p D{(t)- - (t)+)dy = -D{mW) = -DM. (7) 
-b J-b 

Then the instantaneous power P exerted by the plate in self-propulsion, the kinetic energy E of the 
fluid, and the thrust T being the — x-component of the lift acting over the waving plate are given by 

which signify that P is opposite to the work done by specific lift C in making displacement h at the 
rate of /it, E is opposite to the work done by C at the rate of the ^-component velocity W . From these 
integral expressions defining P, E, T, it is obvious that they are related by 

P = E + TU (swimming energy conservation), (9) 

automatically asserting the instantaneous conservation of energy in aquatic animal self-locomotion. 

For the ribbon plate of constant width, the added mass m = pi^lP' is a constant. The integrals 
in dS]) can each be converted to give a sum of two terms, one being a time derivative and the other 
integrable along Q < x < i. For E, it is clear that 

E = m WDWdx='^ {dt + Ud^)W'^dx, E='^u{W\ , (10) 

where the time derivative of a quantity, presumed to be periodic, has the mean value zero (over a long 
time), yielding the integrated value at the trailing edge the time mean of E, denoted by E. 

Similarly, it is straightforward to attain for power P its mean as 



P = mU [htW]^^^. (11) 
Therefore, for z = h{x,t), an arbitrary C^{x,t) function, P,E, and T have their time mean values as 



P = mU { — [ — + U—\ \ , E = —U{[ — + U—\ 
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dt \dt ' ^ dx) y " 2 ^ 1 V ' ^ dx) 

where the last equation results from applying ([9]), and these means are evaluated at x = ^. The foregoing 
presentation in this section is abstracted from Lighthill (1960, 1970) and Wu (1961, 1966). 

These basic relations are of significance in showing that when the mean squares of ht and hx 
in steady cruising remain constant near the trailing edge, the energy shed into the wake, the power 
expended and the thrust acquired all remain constant in time. 

2.1.1. Swimming motion parametric in bodily wave velocity 

Of all the general plate motion z = h{x,t) of ([T]), it is of interest to introduce the special case: 

z = h{x-ct) (0 < x < |y| < 6, t > 0, c G i? > 0), (13) 
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to represent the distally progressing wave with wave velocity c, here abstracted as a one-parameter 
family, parametric in c, for a basic study. By p^ . (jj]) and give 

W = {U-c)h,, C = -{c-Ufh,,, (14a) 

and ()12p becomes (with the overhne (•) omitted) 



P = mUc{c-U)hl, E ='^U{c-Ufhl, T = j{c^ - U^)hl, (14b) 



with hx evaluated at the trailing edge x = i. Here, it is of great interest to notice that in this special 
case, both the mean thrust T and mean power P together become positive (or negative) for c > U (or 
c < U). In particular for this family, there is the hydrodynamic propulsive efficiency, 

TU c + U P-E ^ c-U 

which depends only on the velocity ratio U /c as, the sole parameter and is independent of time, after all 
the common factors cancel out. Further, (jlSp implies that high efficiency can be obtained throughout 
the time for < c - [/ < 2c, e.g., r/ ~ 0.9 for U ~ 0.8c; and r] ~ 0.95 for U ~ 0.9c. 

In general, fishes have other parameters, e.g. longitudinal variations in depth b = b{x) and thick- 
ness, using various appended dorsal, ventral and caudal fins, some of which can even be actively con- 
trolled. These new factors jointly result in differing performances involving additional biophysical effects 
such as in nonlinear and interactive ways that have attracted investigations in broad scope. Here, we 
point out that the propulsive criterion in terms of the speed ratio U/c retains the same standard in 
general as in this simple case of (jlSp . which is not basically affected by other generalizations. 

2.2. Bound and Free Vortex Sheets 

Physically, vortex sheets result from integrating the vorticity vector across the boundary layer 
enveloping the body surface moving through a viscous fluid at sufficiently high Reynolds numbers 
(for the layer being thin). In general, a velocity field in a boundary layer, q = {U + u,v,w), {U 
being the velocity of the irrotational flow enveloping the layer), has its vorticity field a; = V x li = 
{wy — Vz, Uz — Wx,Vx — Uy) = (uji , 0^2 , u}^) , with the local z-axis normal to the body surface (as in ([T]) for 
ribbon plate). Then integrating the {u}i,uj2) distributions across the layer will give the layer-integrated 
vorticity vector 7 = (71,72), which is everywhere tangential to the layer. Making the layer thickness 
vanish on reactive theory, the layer becomes a vortex sheet enveloping the body surface, of strength 
'y{x, y), with V • 7 = d'ji/dx + 8^2 /dy = being the layer-integral ofV-a; = V- Vxu = 0. 

For the waving ribbon plate, integrating the vorticity in the boundary layers on the two sides of 
the plate accordingly yields for \y\ < b, 

/dv V 
-dz = v^-v^ = -2Wix,t)-^, (16a) 
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^dz = u+-U- = —{(/>+ - 0_) = -2Wxix,t)Vb^-y^, (16b) 



satisfying V • 7 = 71^. + 72J, = 0. With W{x,t) given by (|14ap . the vortex-lines can be drawn for 
I7I = ci, C2, • • • over the plate, with each line forming a closed loop because of V • 7 = (see Fig. 1). 
2.3. Vortex Shedding and Wake Momentum 

Having found the instantaneous vortex loops on the sheet bounded to the ribbon plate, everywhere 
regular up to the trailing edge, there encounters the process of shedding the vortex sheet into the 
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ForC>U 

(In the body frame) 



Figure 1: The closed vortex loops on the vortex sheet bound to the ribbon plate are for plate wave speed c greater than 
plate swim speed U . They are shed continuously under the Kutta condition to become free vortices drifting with the fluid 
in the wake. The curly arrows indicate the sense and magnitude of 72, the y-component of the vortex-sheet vector. 



wake downstream. By physical argument, such as so clearly expounded by von Karman and Burgers 
(1934), the flow satisfies the Kutta condition that the velocity field be regular (tangential to the body 
surface at the edge) and flow pressure be continuous at the trailing edge. Once so shed, the shed vortex 
distribution, denoted by 7^, becomes a free vortex sheet and has the properties comprising (i) that it 
satisfies the Helmholtz theorem that it is a material property conserved with the local fluid, i.e. 

Dwlw = {dt + = (or D^c^dt + Udx, if linearized), (17) 

where is the algebraic mean of the velocity above and below the fluid point to which the specific 
wake vortex 7^, will forever be attached. And also (ii), it relates to the circulation T{x,t) defined by 



r(x,t) = <X) V(/)(x,0,z,t) • da; = / 72 (x, t)(ix = — 0_ (0 < x < Xm), 
Jc Jo 



:i8a) 



where the contour C wraps around the plate from the point = {x,0,z = 0_) below the sheet 
clockwise around the leading edge of the plate to point = (x,0, z = 0+) above the sheet, so that 
0+ — 0_ may result with V0 ■ dx = d<j), or ■ dx = (n_|_ — u-^)dx = 72dx, while Xm (> 1) is the location 
arrived at by the starting vortex at time t > 0. Then we have (ii), by Kelvin's theorem, 

/•X PX 

T{x,t)= 72(x,0,0,t)dx = r(x,0) = / 72(x,0,0,0)(ix = (x > x^), (18b) 
Jo Jo 

stating that T{x,t) = for x > Xm since = 0_ there (with no vortex sheets), if T is initially zero. 
This principle is described here for the mid-span of the ribbon plate to also include the 2D case. 

For the ribbon plate waving from at rest, vortex shedding occurs at time t = 0+, since solution 
@ holds at once over the entire plate, with a y-component velocity jump of 71 = V- — v^. From 
(i)-(ii) it follows that the vorticity shedding is continuous at the trailing edge so that behind the trailing 
edge there is a new stretch of free vortex sheet with the same jump in (v- — v+), because otherwise, 
any discontinuity in (v- — u+) would give rise to a singular force impossible to exist. Being free, the 
sheet drifts, by ()17p . with the local fluid of the same velocity, hence no force acting on it, implying 
the pressure continuous across the sheet (in sharp contrast to the bound vortex sheet on the plate in 
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bearing force). More specifically, as the wake gains a circulation 6T as the body moves forward a short 
distance —U5t in the absolute frame, the plate loses, by ()18bp . a circulation by —6T. Furthermore, as 
the lateral momentum mW being shed with a gradient mW^, the free vortex sheet is rotated about the 
y-axis at the rate of dW/dx, with the angular velocity —dW/dx, together with the z-momentum mW 
also rolled into the core of the rolling vortex. This causes a loss of the x-component of momentum by 
—mWWx = —{mW'^/2)x = —dE/dx, which accounts for the energy conservation in the wake. 

For the waving ribbon plate, the wake vorticity 72 (proportional to —W^ by ()16bp ) reaches its 
maximum when the tail (trailing edge) reaches its maximum displacement, lending the free vortex 
sheet to roll at its maximum rate to form the center of a free rolling vortex in counter-clockwise sense. 
Similarly, a rolling vortex forms when the tail reaches its minimum displacement, producing a vortex 
rolling in the clockwise sense. The two rows of trailing vortices clearly pumps the fluid in between to 
gain an excess of momentum going downstream (along a bit undulatory path), thus affording a forward 
thrust needed to overcome the drag acting on the body by the viscous fluid (cf. Fig. 2). This vortex 
street may be called the locomotive vortex street, in sharp contrast to the famed von Karman 'vortex 
boulevard' behind a blunt body being dragged by the viscous fluid. This picture of the wake vortex 
assembled with the waving-pate then accounts for the conservation of momentum in its entity. 



Figure 2: The locomotive vortex street for c > U in the laboratory frame, sending excess momentum downstream and 
giving thrust to the plate (cf. Karman & Burgers 1934, Fig. 105; Lighthill 1969, Fig.8). All the arrows reverse in sense for 
c < U, and drag on the body. 

2.4. Generalizations. 

Aside from the simple case expounded so far, it requires various generalizations to cover most of 
the remaining main factors involved in fish locomotion. Some principal factors raise such issues as the 
body thickness, variations in body depth, db{x)/dx, vortex shedding from dorsal, ventral, and caudal 
fins, with effects on body movement, and vice versa. 

In addition, of the various species of fishes, there are distinct modes of swimming movements 
encompassed by fluid mechanical studies on e.g. the carangiform and thunniform on one hand, and 
those in balistiform and gymnotiform on the other. In the carangiform mode, body undulation becomes 
increasingly more amplified only in the posterior half, or in an even shorter distal end. The mode of 
balistiform and gymnotiform describes a variety of species of fishes that propel themselves by passing a 
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wave along an elongated dorsal and ventral fins with the body held nearly straight. This latter type of 
fish, as told, have been observed to dash out swiftly from a straight tunnel, grabbing their preys, and 
then back straight (without turning) into the tunnel (not wide enough for body to wave), using only 
their elongated fins. Typically, all fish bodies may be regarded as slender. 




Figure 3: (a): Two fish forms show dorsal and lateral fins of mackerel {Scomber scomber) with a row of caudal finlets 
with an operable membrane, (b): A view of balistiform and gymnotiform locomotion by propagating waves along the 
elongated dorsal fin with the body stretched straight. 

2.4.1. Body thickness effect. In general, fish body has a symmetry about its center-plane, which 
may be called the focal plane, across which the body is thin. The irrotational flow past such a body in 
undulation of small amplitude can be evaluated with a distribution of flow singularities (mass source- 
sinks and dipoles) at the focal plane for construction of accurate, even exact solutions by applying a 
scheme known as the focal-plane singularity method, developed for general case by Wu & Chwang (1974), 
Miloh (1974) and otherfl This scheme has been quite actively applied in the naval hydrodynamics 
profession for studies on ship maneuvering and see-keeping in ocean waves, even in designs of bulbous 
bow extending forward from the keel for reducing the self-induced water wave resistance. 

In fish propulsion, the body thickness is found to be secondary in effect by Wu & Newman (1972), 
which can be neglected, or amended if desired, as shown by Newman (1973) and Newman & Wu (1973) 
who employed conformal mapping to map the cross-sectional body contour into a plate for drawing the 
conclusion (see also Yates, 1977; 1983). This also agrees with Lighthill (1975) in giving the added mass 
as m = fiptrlP' parametric in /3, which varies little from unity. 

However, accounting for body thickness is useful to analyzing the effects due to vortex shedding 
by undulatory dorsal and ventral fins such as in studies on the balistiform and gymnotiform modes by 
Wu (1984) and Lighthill & Blake (1990a, b). In these cases, the bodily motion will vary with both x 
and y, as to be discussed in §2.5 and §2.6. 

z = h{x,y,t) {x,y e Sb\ < X < i, \y\<b{x)). (19) 

2.4.2. Vortex shedding from appended fins. This is an important feature on vortex shedding 
from appended dorsal, ventral and caudal fins in fish swimming. For simplicity, we consider the simpler 

^This focal-plane image method can be exemplified by the irrotational fiow of a uniform stream past a sphere represented 
by a certain distribution of a source-sink pair on the front and rear half of the sphere, or a similar pair of stronger magnitude 
on a smaller confocal sphere, or ultimately by a dipole at the focal point at the sphere center. 
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case of locomotion by planar undulation, i.e. z = h{x,t), independent of y as in ([T]) and with body 
thickness still neglected here, but with side fins appended to body of varying depth. This generalization 
can be implemented by dynamical arguments for the anterior, posterior and caudal sections. 

(A) The anterior leading-edge section (0 < x < Xe, db/dx = b'{x) > 0) 

In the leading portion anterior to the widest body section (at x = Xg) of fish body of increasing 
depth with no fin, the lateral swaying movement is analogous to the ribbon plate. It pushes the water, 
giving the transverse slice of water a z-component momentum per unit length in x as in ([6]), 

M=m{x)W{x,t) W{x,t) = Dh = dh/dt + Udh/dx {\y\<b{x)). (20) 

where m{x) = pirb'^ is the added mass per unit length, now a variable. In analogy with d?]), the 
instantaneous specific lift, C, equal and opposite to DA4, is generalized to 

C{x,t) = -DM = -{^ + U-^)[m{x)W{x,t)] {0 < x < Xe). (21) 
ot ox 

This can be used to deduce, in analogy with ([8]), the power, thrust, and energy for this section. 

Furthermore, the total pressure acting on this well-rounded forward body can produce a thrust, 
which has been evaluated by Wu (1961) (1971a, d) along a simple approach treating the body as a focal 
plane with a dipole distribution for the slender-body scheme. The thrust produced due to the cross flow 
across the body section at station x (0 < x < Xe) consists of two parts, Tp + Tg, with Tp given by the 
regular pressure on the plate and Tg the thrust due to the singular pressure at the sharp leading edges 
of the focal plane, as given by Wu (1971a) on plane airfoil theory. 

{'B)Trailing side-edge section (xg < x < Xc, db/dx = b'{x) < 0) 

Posterior to the body section of maximum span (x > Xe), various species of fish (such as dolphin 
fish Coryphaena hippurus, yellowtail Seriola quinqueradiata, mackerel Pleurogrammus azonus, mullet 
Mugil cephalus, etc.) have dorsal and ventral fin edges slanted rearward as trailing edges from which 
the flow leaves smoothly, under the Kutta condition, to form vortex sheets trailing along the body side, 
up to the caudal peduncle at x = Xc from which the caudal fin follows. Under this condition, the effect 
of shedding vortex sheets along body side was found by Wu (1971b, c,d) to make the material rate of 
change in fluid momentum independent of the longitudinal variation of the cross-flow added mass (due 
to body shape alone) in calculating the lift, thus yielding the specific lift in this section as 

C{x,t) = -m{x) DW{x,t) (xe < X < Xc), (22) 

where m(x) = /37r6^(x). This result shows that the effect of vortex sheets shed from side-fin trailing 
edges alongside a body section is to have the momentum transported as if m(x) is a constant, although 
m'(x) < 0. Physically, this implies that the vortex sheets fill in the gap left by the reduced body plate 
width, making the new flow boundary like a ribbon plate uniform in span (of width be = b{xe) V x > Xe). 
We also note that there is no discontinuity between the lift C{xe — 0, t) of (j2ip and the lift C{xe + 0, t) 
of (f22l) provided 6(x) and W{x,t) are both C^(x = Xe) and b'{xe) = 0. 

Lighthill (1970) considered another morphological configuration with the main dorsal fin terminated 
abruptly with a straight transverse trailing edge at x = x^ (xe < x^ < Xc), from which a vortex sheet 
is shed to fill the region between the dorsal and caudal fins (such as in species Cyprinnidae and cat-fish 
Siluroidea). Lighthill's theory shows that in addition to the vortex-free added mass m(x) of a section 
Sx at X (xg < X < Xc), there is also the added mass m{x) due to the frozen body segment, giving 

C{x,t) = —D{m{x)W{x,t) + m{x)w{xs,Tf)} (x^ < x < Xc), (23) 
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where 'w{xs,Tf), being stemmed from (jl7p for material conservation of free vortices, denotes the value 
of the z-component velocity w depending on the fin's trailing edge at Xs and on the Retarded time'' 
Tf = t — {x — Xs)/U, which is the time when the water slice was leaving the fin's trailing edge at Xg and 
reaching the section Sx (located at x) at time t. 

(C) The caudal-fin section (xc < x < £) 

In this section, the interaction between the reentrant vortex sheet and the caudal fin has been argued 
as being capable of augmenting handsomely the lift and thrust generation, especially when they become 
opposite in phase. The problem of evaluating this dynamic interaction is however complicated. An 
earlier attempt by Wu (1971b) contains a deficiency which was removed by Wu &: Newman (1972) who 
took, in addition, the body thickness effect into account, yielding a solution to the caudal fin problem 
in terms of the Abel integral equation which can be evaluated numerically. A consistent slender-body 
theory has been subsequently developed by Newman & Wu (1973), Newman (1973), Newman & Wu 
(1975) by taking into account both the kinematic and dynamic interactions and the effect of body 
thickness on trailing vortices. This problem will be discussed in §2.5 for the generalized case of a thin 
body-fin system moving with displacement h{x, y, t) varying with both x and y. 

To this end, we remark that for fishes known for their high performance using the so-called 'lunate', 
or crescent-moon-shaped tails for propulsion, it calls for abandoning slender caudal fin analysis in favor 
of adopting high aspect-ratio oscillating wing theory. Earlier advances can be referred to Karman and 
Burgers (1934) who showed that high propulsive efficiency cannot be achieved from either the pitching 
or the heaving mode alone, but with the two modes superposed with an optimal phase lag. 

In this respect, it is ingenuous of Lighthill (1970) to introduce a single parameter he called the 
^proportional feathering'' (defined as the product of the maximum incidence angle of a flapping wing 
and the ratio of the flight velocity to wing's maximum tranverse velocity), or -i? = Ua/hu in the same 
notation as above. For heaving and pitching modes, a = kh, k being the wave number, hence i? = U/c, 
implying, by (fTSll . that 1 — <C 1 should be ideal for feathering. Physically, this illustrates how the 
thrust making and its rate of working in the total power expenditure can vary in a competing trend 
after putting geometry, physical principles and simple yet powerful mathematical methods together in 
formulating the control problem to attain valuable results and draw lucid conclusions. It is with this 
clear feathering concept that the 2-D lunate tail wing theory was first expounded by Lighthill (1970). 
This feathering parameter has been effectively used by Chopra (1976) to extend the lunate-tail theory 
to oscillations of large amplitude (for wings of infinite aspect-ratio) and by Chopra & Kambe (1977) to 
predict the hydrodynamical behavior of a variety of lunate-tail planeform of large aspect-ratio. 

2.4.3. Sudden start and change in direction of fish locomotion 

We note that the line of approach just elucidated is also applicable to fish motions with a sudden 
acceleration or change in direction, as found by Wu (1962) for a waving flexible plate swimming with a 
rectilinear arbitrary speed U{t). For acceleration starting from rest, the extraneous power P is found 
positive definite at order 0(t), whereas the thrust can arise only from O(t^); after that the optimum 
motion form can be determined for maximum thrust with power fixed. The problem of fish turning was 
studied by Weihs (1972) for a gold fish turning 90° in a distance of half its length with practically no 
loss in kinetic energy, a feat much beyond any small perturbation scheme. This highly efficient motion 
exhibits itself with the evidence that the locus of the fish centroid appears to be precisely the envelope 
of all the images of the fish centerline (or backbone). Similarly, the problem of fish starting from rest 
was analyzed by Weihs (1973) with the fish centerline set in a preparatory gentle C-shape to send a 
straightening whip of the distal body end to send the fore body dashing forward without turning. These 
two works were discussed extensively in depth by Lighthill (1973a). 
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More recently, this fast-force production problem has been further studied analytically and numer- 
ically using a tail-flapping model by Hu et al.(2004). Both the C-start for avoiding predators and the 
5-start for catching preys rely on placing the backbone in a preparatory C oi S shape for shedding an 
impulsive tail vortex as found to be important to impulsive thrust generation. The efficiency can be 
optimized by taking larger deflections with no delay at stroke reversal. An integrated method has been 
developed by Yang et al. (2008) to study 2-D self-propelling aquatic and aerial locomotion involving 
unsteady body deformation dynamics, obtaining results in close agreement with experiments. 

2.5. A generalized slender-body theory of fish locomotion 

The general case of body- fin movements that vary with both x and y, as that found in balistiform 
and gymnotiform, has been investigated by Wu (1984) and also by Lighthill & Blake (1990a, b). This can 
be taken as a boundary- value problem based on the linear slender-body theory, with the flow represented 
by a vortex distribution on the body focal plane or using conformal mappings. This irrotational flow in 
a fluid incompressible and inviscid has the velocity potential (f>o {x,y,z,t) such that 

^, = Ux + <l){x,y,z,t) (v2,/.o = = 0), (24) 

giving the perturbation velocity (relative to the mean body frame of reference) as 

u = iu,v,w) = V ^ {V = {d,,dy,d,)). (25) 

For slender bodies of length i and depth b characterized by slenderness 6 = b/i <^ 1, we can scale 
y, z hy b and x by £, thus converting the equation = in (f24|) under the new scales to 

(t>yy + (t>zz + S'^(t>xx = 0. (26) 

If the body motion displacement is also of order z = h{x, y, t) = 0{6i), ()26p then provides the foundation 
for developing the slender-body theory starting with the expansion (j) = (j)Q + 5'^<j)i which produces, 

by substitution in (|26]l . a set of ordered equations for 0^ (A; = 0, 1, 2, • • • ) as 

(pyy + (pzz = {k = 0), (t>{k)yy + <t>{k)zz = 4>(k-l)xx {k = l,2,---). (27) 

In practice, the leading order may suffice unless h is of large amplitude, as we will discuss in §2.6. 
In the leading order, the Bernoulli equation for fluid pressure p can be linearized to give 

D0 = $, <^ = {p^-p)/p {D = d/dt + Ud/dx), (28) 

where ^> is Prandtl's acceleration potential, p the fluid density. From (27) and (28) it follows that <I> 
is also harmonic, with its conjugate function ^. In terms of complex variables as used by Wu (1971b, 
1984), the coordinate C, velocity potential /, acceleration potential F, and complex velocity X) 

( = y + iz, f = (f)^i^^ F = ^ + i'i!, x = v-iw = df/dC,, (29) 

are analytic functions of one another and are related by 

Df = F Dx = dF/dC {D = d/dt + Ud/dx). (30) 

The boundary conditions for this problem can be expressed in terms of either or $ as follows. 
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I. Boundary conditions in terms of (p II. Boundary conditions in terms of $ 

{i) {(t>,)± = Dh = W{y;x,t), = DW {y;x,t) ((x,y)G5fe), 

{ii) = 0, $± = 0, {{x,y)£S^), 

(Hi) cf>± = 0, ^'± = 0, {{x,y)eSc), (31) 

(iv) = 0, $± = 0, ((x,y) G T.E.), 

(v) f = O(r'), F = 0{C'), X = 0{C') (ICI^oo). 




x^o a x^x, b x^x, c X-/ 




Figure 4: The various sections of flow regions for analyzing fish propulsion: (A) anterior leading-edge section (0 < a; < Xe), 
(B) trailing side-edge section {xe < x < Xc), and (C) caudal fin section {xc < x < £). 



Here, the kinematic conditions I(i) and II(i) over the body planeform Sj, at z = prescribe the fish 
body lateral sway with velocity W{y;x,t) (V \y\ < b(x)). In (ii), Su, is the wake planeform resulting 
from projecting the vortex sheet shed from trailing edges of side fins onto the z = plane. In (iii), Sc 
denotes the region at the z = plane which is complementary to 5^-1-5'^; this condition is due to the 
symmetry in (i). The Kutta condition (iv) requires the pressure to remain continuous at the trailing 
edge (T.E.) of side fins and caudal fin (with the pressure p gauged by pbj) as in (ii) and (iii)). And the 
asymptotic behavior in (v) stems from Kelvin's circulation theorem so that the flow solution be source 
free and vortex free at = oo. Schematically, this is shown in Fig. 4. 

The problem formulated above is in a form of the Riemann-Hilbert problem. By classical analysis, 
the required solution can be obtained for different longitudinal sections 5^^ of the body as follows. 

(A) The anterior leading-edge section (0 < x < x^, db/dx = b'{x) > 0) 

In this section, the side edges along y = zizb{x) are assumed well rounded and moving well feathered 
so that the velocity field of the plane flow at station Sx can be represented by a distribution of vorticity 
7(x,?/,t) on line L at body section —b<y<batz = 0, i.e. 

X = v-iw = ^ [ ^^^^^dy^ (C = y + iz^Sb). (32) 

For condition I{i), we let C ^ ± ^0 {\y\ < b) on Sb and apply Plemelj's formula, giving 

v±- iw± = ±\-i{x,y,t) + ^.P I ^^^'^^'^^ dyi {y,yi e Sb), (33) 
2 2Tn yi-y 

where the integral designated with P assumes Cauchy's principal value. Hence, 'y{x,y,t) = V-^- — V-, 
and by I{i), we have the integral equation for "f{x,y,t) given as 

W{x,y,t)=w± = ^P [ 2lMlll)dyi^G[7(x,y,t)] {\y\<b), (34a) 
2vr Jl yi-y 
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where G denotes the integral operator as shown, and we have the solution (P for J being omitted) as 



2r H(yJI_Wi^ 

71" Jl H{y', h) y' - y 



where H{y, h) = \j\P^—]p' is the complementary solution with the symmetry (i) of (|3ip . and is the 
inverse operator to G 9 G~^G = GG~^ = 1. This is readily verified, for substituting ()34bp in (|34ap 
results in an identity for W{x,y,t) by virtue of Poincare-Bertrand's formula: 

[ dyi f f{yi,y') ^ , 2 \ , [ a ' [ fivi^v') ^ (oa \ 

/ / — ; dy = -TT f{y,y)+ / dy / — -dyi. (34c) 

jLyi-yJi y -yi Jl Jl [v -yi){yi-y) 

Therefore, the velocity field in this section results from (|32p and ()34bp as 

% = xiC) = v-iw = ^ [ -^G~^[Wix,y,t)] (0<x<xe). (35) 
dC 27r? Jl yi - ( 

(B) Trailing side-edge section (xg < x < Xc, db/dx = h'{x) < 0) 

Within this section, the slant elongated dorsal and ventral fins are trailing edges (with b'{x) < 0), 
and with the fins waving in displacements that may vary also with y, the Kutta condition being a 
prerequisite. Under this condition, we note that there is a complete analogy between the two sets of 
boundary conditions in I(i) on ^2 = Re{df /dC,) = W and II(i) on = Re{dF/dC,) = DW, other 
conditions being equal. Whence, by analogy with (j3ip . we attain for dF/dC, the solution: 



dF _ I f dyi 



[ -^^G-^[DW{x,y,t)] {xe<x<x,), (36a) 

Jl yi-C 

where operator G~^ is given in ()34bp . Upon integration, we obtain F under condition (v) as 



f\ 1 /"S^ r ^2/1 f' H{yi) DW{x,y',t) ^ , 

F{C,x,t) = ^ dC -/ dy (36b 

^^^^ J-oo J-b yi-C J-b H{y') y' - yi 

which is seen to satisfy the Kutta condition in particular, since F is continuous everywhere, including at 
fin's trailing edges at y = ib6(x), by (iii) in ()3ip . This F field provides the pressure distribution within 
this trailing-side-edge section. The corresponding /(^) can be obtained, straightforwardly by integrating 
(30), i.e. Df = ft + U fx = F, along the linearized characteristic lines x — Ut = ^ = constant, starting 
from X = X(, (at which / is known from (j3ip ). with F given by ()36bp . Then the velocity field can be 
deduced from x = df /d(, thus yielding the (±)-side limit of x = i' — iui over Sb and Sw (for the details, 
see Wu, 1984). The result of w± = w{y, ±0; x, t) = Wy{y; x, t) gives the \sidewash' induced by the entire 
vortex system consisting of the bound vorticity lying in the planar body surface Sb and the free vortex 
sheet over S^, : {b{x) < \y\ < be, Xg < x < Xc). This induced velocity field, when convected with the 
free vortex sheet into the caudal fin section, is responsible to assuming an important role in interacting 
with the caudal fin to enhance its thrust production and propulsive efficiency (cf. Wu 1984). 

The merits of using the acceleration potential $ are clearly exhibited here in the analysis showing 
that <I> vanishes on all wake S^, and complementary Sc regions including, by (iv) of ()3ip . the trailing 
edges. Only when the velocity field is finally attained by integration along the characteristic lines do 
the velocity jumps appear across the free vortex sheets on and with velocity being continuous on Sc- 

(C) The caudal-fin section (xc < x < £) 
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In this section the caudal fin has new leading edges in presence of the vortex sheets incident from 
upstream which in turn interact with the fin motion. In the framework of linear theory, the velocity 
field can be constructed by superposition of two components, one being Xv which is that induced by the 
entire vortex sheets and is convected intact downstream (with velocity U) into the caudal fin section, 
and the other, Xc: which is the new velocity due to the motion of the caudal fin itself, giving 



x(C; X, t) = xviC, X, t) + xciC; x, t) {xc<x < e), 



(37) 



where x = i marks the caudal-fin trailing edge. Based on the known sidewash velocity, Wv from (B) 
above, induced by the vortex sheets on the caudal-fin surface (which may be slender or a lunate tail of 
large aspect ratio), we can obtain for Xc its solution in integral form (cf. Wu, 1984). 

Here, we note that if Wy vanishes and if the caudal fin is slender, this Xc reduces to (j35p for the 
anterior section. When ^ 0, the interaction between the caudal-fin movement and the vortex sheet 
sidewash is seen possible to enhance greatly the thrust generation and propulsive efficiency, especially 
when their own induced velocity fields maintain opposite in phase because the two sidewash terms 
pertaining to Xv and Xc then totally augment each other. Typical cases of numerical examples have 
been investigated by Su et al. (1983) and by Yates (1983), yielding results that lend quantitative support 
to the qualitative prediction. We refer further details also to Wu (1984, 2001b). Finally, we remark 
that if W is independent of y, then (j32]) reduces to the simpler solution as first expounded in §2.4. 



2.6. Large-amplitude Fish Swimming 

For the important problem that fish undulatory motion can be large in amplitude, Lighthill (1971, 
1975) pursued this generalization with two principles, (i) Kinematically, the trajectory of a fish moving 
with arbitrary amplitude can be described by using one Lagrangian coordinate, denoted here by ^, 
initially stretched straight along the fish backbone jointly with the local cartesian system so that at 
time t a point x = (x, y, z) = x{^, t) on it is prescribed relative to an absolute frame, with fluid at rest 
far away, in which a;(^,0) = ^ (0 < ^ < 1) at t = 0. The fish is assumed to move in the horizontal 
y = plane tangentially along a trajectory (x(^, t), 2(^, t), t), with unit tangent r = {dx/d^,dz/d^) 
and normal n = [—dz/d^, dx/dS,), both lying in the xz-plane with the fish depth being always vertical. 
The backbone is further assumed inextensible, i.e. (dx/d^)'^ + (dz/d^)'^ = 1. (ii) Dynamically, it travels 
forward with cartesian velocity u = {dx/dt,dz/dt), or with tangential velocity u forward and sidewise 
with normal velocity w, 



dx dx dx dx ^^.^ , 



w 



u ■ n 



dz dx 
'di'di 



dx dz 



WiC,t). 



(38) 



This normal velocity W gives at the local ^-station the lateral momentum per unit length equal to 
mW (m ~ /97r6^(^)), b being half the local fish-depth. With this mW, it is then possible to give the 
instantaneous reactive force F = {P,Q) exerted by fluid on fish, after Lighthill (1975), as 



F = {P,Q) 



dz dx. 1 ^^^n 

-1^7' +1^ + -mW^n 
dt dt 2 



(39) 



with [S^ = £] indicating the value given at the caudal end. Here, —F is the force exerted by the fish on 
the fluid in propulsion. The various terms in ([39|) generalize those attained in linear theory. 

This important theory can be further developed here to cover the general case with vortex shedding 
from appended fins moving at large amplitude also taken into account. The underlying principle is to 
take a local instantaneous cartesian coordinate system x = (x, y, z) fixed at the fish centroid at time t, 
with its centroid moving forward tangentially to its backbone. Then the local solution for the velocity 
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potential (/'(C(C> 0> ^-nd acceleration potential ^{Ci(,, t),t) can be derived for the anteroior and posterior 
sections by complete analogy with the solution (j35p for (p + iifj and with (|36p for F = ^ + so that 
the free vorticity distribution can be all evaluated, and the final solution follows by time marching. 

It is gratifying that Lighthill's classical works on fish locomotion have been nicely extended by 
Singh and Pedley (2008) for analyzing finite fish movements and by Cheng (1998) in developing the 3-D 
flow models to implement 'recoil' corrections in order to ensure zero net force and moment acting on fish 
body. Jointly, a 3-D boundary-layer algorithm has also been developed using a vortex lattice method 
with applications to various cases for comparisons with analytical results available. A significant step 
has been further ventured by Pedley and Hill (1999) in integrating fluid mechanics and fish internal 
mechanics for insight investigations of large-amplitude undulatory fish swimming. 

3. BIRD FLIGHT 

Theory of oscillating airfoils has important applications, including that applied to studies on bird 
and insect flight and perhaps more in practice to aero-elasticity for expounding the fluttering forces and 
effects on airplane wing- fuselage system flying through air with turbulent patches and sharp gusts. 

Of the various linear theories developed in the 1920s, Wagner (1925), in particular, first derived 
an integral equation of great significance to provide accurate determination of the starting vortex shed 
from an airfoil moving smoothly or impulsively forward. For a 2-D airfoil z = h{x,t) ( — 1 < x < 1) in 
oscillation, Wagner's integral equation reads 



where ro(t) is the quasi-steady value of the circulation of the wakeless airfoil determined on steady 
airfoil theory with time t frozen as a parameter, and jw{x,t) is the free vorticity in the wake (x > 1) 
with the starting vortex arriving at Xm{t) for t > 0. Otherwise, these linear theories are thought to 
be either very complicated or hindered from being generalized, and not very transparent in physical 
conception. These shortcomings have been overcome by von Karman and Sears (1938) who created a 
simple and elegant restructuring of the bound and free vortex sheets varying on Kelvin's principle of 
conservation of circulation, or the moment of momentum, as will be further elucidated in §3.3. 

3.1. Oscillating rigid airfoil. 

During the early stage, rigid airfoils can oscillate in heaving and pitching as the only two mode^. 
In this respect, Karman and Burgers (1934) noted that high propulsive efficiency can be achieved only 
with the two modes superposed with an optimal phase lag as the single parameter. Various attempts 
followed to combine the two monochromatic modes with differing phase lag, such as that of seeking one 
phase lag for the highest propulsive efficiency with thrust specified. 

Alternatively, Lighthill (1970) proposed to take the pitching axis location as a new parameter, but 
with the heaving phase-lagged by 90° behind pitching, i.e. 



where h is the heaving amplitude and a the amplitude of pitching about the axis at x = 6 as the new 
parameter. With (41), the energy expelled to the wake was found to have a sharp minimum for b = 1/2 
(at the 3/4-chord point) whereas the rate of working increases with increasing b, hence the optimal 

■^Separately, heaving alone has its propulsive efficiency rjh falling with increasing reduced frequency a = toc/U, {lo being 
the circular frequency, c half the wing chord) from 77^ = 1 at cr = to its asymptote at rih ~ 0.5, whereas by pitching 
(about the mid-chord), its efficiency rjp rises above zero at a = CTo = 1.781 and increases to rip ~ 0.4 (Wu 1968). 




(40) 





(41) 
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position is slightly behind the x = 3/4-chord point (for details, see Lighthill, 1970, p. 297). This was 
in one of the manuscripts mutually communicated then between Lighthill and Wu for discussion and 
reference. So (41) was used early by Wu (1971b) to study the tail movement of a porpoise in cruising 
with data reported by Lang and Daybell (1963). The tail fluke angles relative to the path of tail base 
predicted by Wu (1971b) and the measured are found in excellent agreement with the pitching axis 
found at a mean of the 0.793-chord point, well in the favorable range predicted by Lighthill. 

3.2. Generalizations. 

One important new issue is to consider wing flexibility for unsteady changes in its camber shape by 
command. To explore the resulting effects on propulsive qualities, locomotion of a flexible wing was first 
studied by Wu (1961). For an airfoil starting with acceleration from at rest, the unsteady effects have 
been examined by Wu (1962); the small time solution exhibits the decay of transient terms, enabling 
the optimum shape determined for the maximum thrust with fixed power. 

There are several salient features of significance in the world of self-locomotion of aquatic and aerial 
animals using lifting surfaces such as wings and body-appended fins. First, the wings are in general large 
in aspect-ratio, a feature that would suit for an unsteady lifting-line approach. Secondly, the periodic 
flapping of the wing generally involves changes in surface-shape function, e.g. from a stretched-straight 
wing in downward pronation stroke to a form with a noticeable arched camber and spanwise bending in 
upward supination stroke. Further, in swift maneuvering, the wings may bend and twist asymmetrically 
to change in orientation and trajectory, e.g. in the beautiful performance of a humming bird using a 
figure-eight wing flapping in hovering its body flxed in front of a flower, and then suddenly fleeting off 
in a flash. All these features are so strongly nonlinear and time-dependent that a comprehensively valid 
theory would have to take all these factors fully into account. 

3.3. A Fully Nonlinear Unsteady Theory for a Flexible Wing 

Recently, a fully nonlinear and unsteady theory has been developed by Wu (2001b)- (2007) for 
evaluation of unsteady flow generated by a 2D flexible wing moving in arbitrary manner for modeling 
bird flight and flsh swimming. This is illustrated here for facilitating applications. 

Kincmatically, the wing shape function and the resulting wake formation consists of a 2D flexible 
wing Sb{t) of negligible thickness, moving with time t >0 through the fluid in arbitrary manner, as can 
be specifled parametrically using a Lagrangian coordinate system (^, rj) jointly with a local cartesian 
system x = (x, y, z) for the movement of a flexible wing with arbitrary variations in wing shape and 
trajectory to produce an irrotational flow fleld in an incompressible and inviscid fluid. This identifles 
a point X = X{^,t), y = Y{^^t) on the boundary surface S{t) = Si,{t) + Sw{t) comprising the body 
surface >S'(, and a wake surface 5^, with S{t) lying at time t = over a stretch of the ^-axis (at 77 = 0), 
X{^, 0) = ^ (—1 < ^ < 1), y(^, 0) = 0, and moving with time t > in terms of the complex coordinate 
z = X + iy = Z{^, t) as 

Z{C,t) = X{^,t)+iY{^,t) on Sb{t) : {-1< ^ < 1) + S^t) : {1< ^ < U), (42) 

with ^ = — 1 at the leading and ^ = 1 at the trailing edge of the wing. From the trailing edge a vortex 
sheet is being shed smoothly to form a prolonging wake Syj, and identifies the path Z{S^rn,t) of the 
starting vortex shed at t = to reach = ^mii) at time t > 0, with the fluid initially at rest in an 
absolute inertial frame of reference. The flexible Si,{t) is assumed inextensible {\Z^\ = \dZ/d^\ = 1, or 
X| -I- = 1, 1^1 < 1) and the point ^ on Si,{t) moves with a prescribed velocity W{^, t) = U — iV, 

W{(, t) = U -iV = dZ/dt = Xt- iYt {\^\ < 1, t > 0; Z = X - iY), (43a) 
WdZ/d^ = {X^Xt + Y^Yt) - i{X^Yt - Y^Xt) = U, - iUn, (43b) 
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Us{(,,t) being the tangential and the normal component of the wing velocity as prescribed, and 

similarly for the wake surface 5'u,(t) for (1 < ^ < £,m), but with its velocity derived with the solution. 




Figure 5: (a): The Lagrangian coordinates (^,77) describe arbitrary motion of a 2-D flexible wing moving along an 
arbitrary trajectory 77 = through unbounded fluid at rest at infinity, (b): The wing moves with (i) translational velocity 
(dZo/dt) at incidence angle a{t), (ii) rotational angular velocity ^l{t), and (iii) unsteady camber function Z{^,t). 

More specifically, the wing Sb moves with a shape function Z{^,t) consisting of (i) rectilinear 
translation with velocity dZ^jdi at incidence angle a(i), (ii) rotation with angular velocity and 
(iii) flexing with camber function Z(^, t) which is prescribed by (see Fig. 5, the inset) 

ZiC, t) = Zo{t) + e'' Z{i, t) Z{i, t) = X{i, t) + iY{i, t) = X(e, t) + iF{X{i, t),t), (44) 

where Z{l,t) — Z{—l,t) gives the chord line inclined at angle 9{t) relative to the x-axis fixed in the 
absolute frame, with pitching at angular velocity n{t) = —dO/dt (> for nose-up pitching). Here, Z{(^,t) 
is the complex camber function given in the body frame, with the origin set at ZQ{t), corresponding to 
X{0,t) = 0, Zq = Z{0,t) — ie'^y(0,t), which is the projection of point Z{0,t) onto the X{^,t)-axis, 
with Y{^,t) = F(X{$^,t),t) being the real camber function, and with \Z{l,t) — Z{—l,t)\ = a — (—6) = 
{a + b) = c{t) the chord length. 

For the best resolution, the simple and elegant physical concept crystallized by von Karman and 
Sears (1938) has provided an ingenious restructuring of the vorticity distribution over S{t) as follows: 

onSbit): 7(e,i) = 7o(e,t) + 7i(e,t) (-K ^ < 1), 

on5^(t): litt) = 7witt) {l<C<Cm)- 

Here, 'yo{^,,t) is the vortex bound to in the wakeless "quasi-steady" flow past which can be 
determined by steady airfoil theory for the original prescribed Un{^,t), with time t frozen to serve 
merely as a parameter for the moment. Vorticity 71 ('^,t) is the additional bound vortex induced on Sh 
by the trailing wake vortex 7«,(^,t), yet unknown, such that 71 and 7^ jointly bear no change to Un 
over Sb so as to reinstate the original time-varying normal velocity Un{C,t) prescribed on Sb{t). 

Briefly, a cambered wing sets in motion at t = 0+ in a uniform stream of velocity U produces a 
quasi-steady velocity field represented by a distribution of vorticity 70(^1 0+) at Sb{t = 0+), i.e. 

w{z, 0+) = u-iv = ^ [' 7°^^'°^^ {z = x + iy^ 5,(0+)), (45) 

27^^ J_i Z(^, 0+) - z 

where Z(^, 0+) is an initial cambered foil across which w has a jump given by Plemelj's formula as 

uf - ^ut = «^^(e,0+)^ = ±^7o(e,0+) + I' di' (Z(^,0+) E 5), (46) 
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with velocity Us tangential and Un normal to 5*^(0+) and with the integral assuming Cauchy's principal 
value. Here, ([l6j) gives 7o(C,0+) = uf — uj (—1 < ^ < 1) and = u~ = Un normal to S. With Un 
given on S^,, 70 satisfies, for all {CO £ S^,, the integral equation: 



For a flat wing, the residual kernel giC^t) = due to - Z{^) = e*''(*)(^' - ^) for arbitrary e{t). 

For wings of small camber, solutions for 70 converges rapidly by iteration (cf. Wu 2007). 

At the same time, a vortex sheet element jwiS,,ii) is shed into the wake of a small length 
= UAt lying in 1 < ^ < 1 + A^, at this first step in time-marching calculation with t = 
ti,t2--- , ifc-ii ifc, ■ ■ ■ (with tk — tk-i = At taken aptly small). This vortex sheet Jw{£,,ti) induces 
on Sh a velocity distribution as uf iCh) - iuniiM) = w^{i,ti)dZ/di (1^1 < 1), similar to (06]) for 
t = 0+, of which the normal component ttn(C)*i) (1^1 < 1) has to be canceled by a new bound vor- 
tex 7i(^,ti) on Sh{ti) in order to maintain that Un varies as specified for the plate motion. And 
7i(^)^i) (l?l < 1) can be solved just like 70(C) *i)- Thus, the first step involves only one unknown since 
7i is a function of 7^^. Finally, 7«;(C, ti) can be solved by applying the Kelvin Theorem requiring that for 



t = ti,To + Ti+r^ = (70(^,^1)+ 71 (C,ii))c^C + /i ' 7wiS.,ti)d^ = 0, and similarly for t = t2,h, ■ ■ ■ 



in time marching for the solution. 

Carrying out the time steps to t = (A; = 1, 2, • • • ) then yields an integral equation for as 



This is Eq.(37) in Wu (2007), with the nonlinear terms Ni){(^,t) and A^^(C,t) given therein, and is 
proclaimed as the general wake-vorticity theorem in terms of which nonlinear integral equation for wake 
vorticity 7^ to satisfy exactly for the solution to be exact. It generalizes Wagner's integral equation 
for the linear case to account here for a fiexible wing varying in arbitrary wing shape and its arbitrary 
trajectory. In the linear limit, both N^j = and Nf^ = 0, hence (|48p reduces to Wagner's linear integral 
equation (|40]) . 

3.4. Comparison between theory and experiments. 

This nonlinear theory has been applied by Stredie (2005) and Hou et al. (2007) to carry out 
computations of various unsteady motions of a 2D airfoil to attain results of high accuracy in all the 
cases pursued over a diversified broad scope. Of them two special computational studies are presented 
here for comparisons with two well noted wind-tunnel experiments. The first is for an airfoil performing 
heaving oscillations at reduced frequency a = ojc/U {uj being the circular frequency, c the half-chord, 
U the free stream velocity) for the one-parameter family solution parametric in a and computed for 
a = 2.0 and heaving amplitude h = 0.038 as used in the experiment by Lai (2002). The numerical 
results (for time step At = 0.001 to 2630 steps) is shown in 

Fig. 6, with the central photo depicting the observed wake vortices shed from the trailing edge situated 
near the left border, growing in size with their senses indicated on the top, and with the corresponding 
numerical result shown in equal scale right below with a single vortex line as the centerline of the 
real physically diffusing vortex wake. The qualitative and quantitative agreement between theory and 
experiment is excellent. More specifically, the distance between the centers of the second pair of vortices 
is measured in the photo to be very nearly 0.4 chord, so is exactly the computed result. On the right 
of Fig 6 is a computer plot of the longitudinal fluid velocity, in excess of C/, produced by the heaving 




(47) 




,m — 




(48) 
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Time steps 

Figure 6: (a): Heaving wing; Upper panel: rotational sense of eddies; Mid panel: experiment by Lai (2002); Lower panel: 
Wu's numerical results; (b): The flow velocity of a heaving wing, recorded at 0.4f chord downstream of the trailing edge. 




Figure 7: A wing pitching about the 1/4-chord point and moving from left to right; Upper panel: Wu's numerical results; 
Lower panel: experiment by Koochesfahani (1989). 



wing at a point 0.41 chord downstream from the traihng edge. The negative values (for fluid flowing 
upstream) in the first half period is an interesting recording of the local velocity field during the shedding 
of the very first starting vortex, then followed by periodic fluctuations of the local wake velocity, which 
is largely positive, thus implying a forward thrust exerted by the fluid on the heaving wing. 

The other case is for an airfoil performing pitching oscillations, about the quarter chord (or about 
the mid-cord plus heaving) computed with those data used in the experiment for the pitching amplitude 
of a = 2° with time step At = 0.0005 to 3000 steps. Figure 7 shows the comparison of the numerical 
results with the experiment accomplished by Kooschefahani (1989) in which the photo of the wake 
vortices shed from the wing (moving from left to right) is shown below the corresponding numerical 
result. The comparison is again excellent in precision. Both comparisons stress the utmost importance 
in having the starting vortex computed in the very first time step in high accuracy, for a small artificial 
error is found to lead to growing departure of the wake vortices from the accurate results. This also holds 
for all the cases examined, including mixed heave and pitch, time-periodic bending of the wing plate 
for simulating membrane wing (like bat flight), impulsive start in incidence and camber of a Fourier 
flexible wing, etc., as shown in details by Hou et al. (2007). 
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The nonlinear unsteady motions of a 2D flexible wing has been extended by Hou et al. (2006), to 
3D for a rectangular plate advancing normal to the flow. These 3-D numerical results are found also in 
excellent agreement with the corresponding experiments (Ringuette et al. 2007). The generalization to 
3D flexible wing of arbitrary plane form could adopt some lifting-line approach for aspect-ratio no less 
than four, such as for the case with given spanwise free stream distribution as shown by Tuck (1975) or 
other types of lift-line approaches. 

4. INSECT FLIGHT 

The importance of studies on insect flight is evident. The diversity in all aspects of the insect 
kingdom is indeed overwhelming, with close to 10® species of winged insects, all self-propelling by 
flapping their wings of all sorts of morphological varieties, exhibiting various aerodynamical properties 
and distinct features as reviewed by Triantafyllou et al. (2000) and by Wang (2005). To our mankind, 
one might wonder what would have become to the plant world without honey bees and butterflies, 
and of course there are pestering insects too. For large insects like moths, the two wings may span to 
5cm. For medium to small insects down to 3mm, the wings are thin plane, corrugated, supported and 
executed by a somewhat elastic wing root to a single pair of wings (like in bee) or two pairs of forewings 
and hindwings (like in dragonfly), with the aspect-ratio varying from 3.8 to 12. The wings of very small 
insects such as the various thrips can have deep fringes composed of bristles (setae) projecting outward 
from wing margins to form slotted sheets and to passively propagate waves generated by wing flaps so 
effectively that they can outperform manyfold those achievable by tiny airfoils of equal size flying at 
the same Re, as reported by Kuethe (1975). 

4.1. Principal Features and Fundamental Properties 

Kinematically, insect wing flappings can appear in various modes, e.g. in periodic pronation down- 
strokes and supination upstrokes, or in hovering flight, with body fixed in position, by flapping the pair 
of wings in repeated figure of 8 (or figure oo). The wing pair can either keep flapping apart (without 
contact), or with periodic contact in the clapping mode. For the latter case, the pioneering contribu- 
tion to novel insect flight is the 'clap-fling mechanism' discovered by Weis-Fogh (1973, 1975). This 
mechanism, clap to fold together and fling to fly apart, has been extensively investigated. 

Aerodynamically, there are two basic parameters in the wing flapping operation. One is the 
Reynolds number. Re = Ui/u, U being a typical velocity, £ a typical length, v the kinematic vis- 
cosity of the air, which is found to fall in a range of meso Reynolds number^ like 10^ — 10^. The other 
is the Strouhal number, St = uj£/U, which represents the unsteady effects due to flapping of circular 
frequency w. In addition, an outstanding feature in wing stroke has been observed that the insect wings 
in general are found to sweep at large incidence angle (of attack), like in the 20° — 40° range, so much 
in sharp contrast to bird flight that its mechanism and effects inspire exposition. 

To this end, we may cite a famed study on the locust Schistocerca gregaria in a role-model collab- 
oration by Weis-Fogh and Jensen (1956), a zoologist and an aerodynamicist of renown. Placing in their 
wind-tunnel a specimen of weight 2g, wing span 4cm, chord c = 0.75cm, in a wind of speed U = Am/ s at 
Re = Uc/i^ ^ 2000, they found the wing beat frequency at / = 20Hz (uj = 2i:f) or the Strouhal number 
St = ujc/U ~ 1/4, which is considered very high for the unsteady effects, as can be attested by the 
Stokes boundary-layer thickness 5 = \j2v juj being small (like that for a flat plate oscillating in its own 
plane in air at circular frequency w), which in this case gives b = 0.48mm, a negligibly thin boundary 
layer on the wing. This finding then establishes the first basic principle (i) that for St ^ 1/Re, the air 
can be assumed incompressible and inviscid for aerodynamic evaluation of lift, power and efficiency of 

^The term of 'small (or low) Reynolds numbers' Re had been commonly used, prior to the insect-flight research days, 
in both mechanical and biological literatures on microorganism locomotion with 10~^ < Re < 10~^. 
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insect flight. This assertion should even be enhanced for smaller insects, because in general, the smaller 
the insect, the lower the Reynolds number, but the higher the beating frequency, hence the higher the 
Strouhal number. As a result, the lift coefficient Cl gained by small insects in hovering flight can be as 
large as from 4 to 6 (Weis-Fogh 1973) due to the unsteady high lift by insect wing flapping. 

An additional issue is concerned with the high incidence angles held by sweeping insect wings. This 
actually is related to the clap-fling mechanism, first discussed by Lighthill (1973b), later modified by 
Maxworthy (1981) who discovered flow separation from the leading edge of a model wing in forming 
vortex sheets moving towards the wing tip as found from his laboratory experiments. In this regard, 
a comprehensive view has come to focus as follows. Ellington et al. (1996) observed flow separation 
around the leading edge of wings of tethered hawkmoth and a flapper model, forming a strong and stable 
leading-edge vortex (LEV) during downstroke to lend high lift forces. Dickinson et al. (1999) measured 
unsteady forces acting on a fruit-fly wing model hovering at Re = 136, / = 14:5Hz, finding out three 
underlying mechanisms: delayed stall in down- and up-strokes, rotational circulation and wake capture 
coming with rapid rotations of the wing in stroke reversal, which can provide 1/3 of the mean lift. 
More recently. Sun & Tang (2002) numerically analyzed the 3D N-S equations for the hovering model, 
finding that flow separates from both the leading and trailing edges under the Kutta condition to form 
the leading-edge vortex plus its image (LEV&IM) in the wing-body system and the trailing-edge vortex 
plus image (TEV&IM) as the three major mechanisms. Thus, it is clear that the high incidences held 
by insect wings are responsible for both (LEV&IM) and (TEV&IM), with the latter keeping the former 
stable in float, as noted by Yu, Tong &; Ma (2003). These advances in conception have stimulated 
interest in insect flight to the highest level of activities in the last decade. 

More recently, developments have been made in aerodynamic theory and modeling, computational 
and experimental studies on various aspects of insect flight for our expository review below. 

4.2. Aerodynamic Theory and Modeling of insect flight 

It has been noted that effective aerodynamics at Re < 10^ is of great technological interest and 
a fundamental scientific challenge. Also noted, by Lian et al. (2003), the aerodynamic performance 
of a wing deteriorates considerably as the Re number decreases from 10^ to 10^. For the unsteady 
leading-edge vortex (LEV) formation and shedding of wake vortices in insect wing flapping, Ansari et 
al. (2006) developed two novel, coupled, nonlinear wake-integral equations for evaluation of force and 
moment for 2D inviscid flows. Separately, Bandyopadhyay (2009) regarded the effects due to the LEV 
phenomenon being responsible for unsteady high lift. 

A fundamental theoretical modeling for wing flapping of small insects has been developed by Yu, 
Tong, &; Ma (2003) for a 2D flat plate oscillating in heave and pitch, shedding vortex sheets from both 
the leading and trailing edges, as well as with swift rotation of the viscoelastic wing roots at stroke 
reversal. Employing flow singularities of dipoles and vortices on the plate to construct the potential flow 
in an incompressible and inviscid fluid, results obtained have simulated the experimental conditions of 
Dickson Sz Gotz (1993) for the unsteady lift coefficient Cl at incidence angle of a = 31.5°, Strouhal 
number St = 0.35, which are found in good agreement. Their is given in three components, one 
for the added mass, the other two for LEV&IM and TEV&IM each, exhibiting, interestingly, that the 
TEV&IM is dominating in general. Also, their results for a pitching wing are found in good agreement 
with the CFD results by Hamdani &; Sun (2000). This work was revisited by Yu et al.(2005) for further 
expounding on the inertia effect, the LEV and TEV. These works have been adopted by others for 
modeling insect flight in applications. 

In a valuable review, Tong &; Lu (2004) discussed several noticeable issues concerning insect flight. 
One issue points out some shortcomings in quasi-steady schemes for unsteady insect flight that appears 
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in need to open a new field for unsteady aerodynamics at relatively high Strouhal numbers. Another 
issue stresses the fact that the pronation down- and supination upstrokes are along inclined down- 
and up-paths rather than the commonly assumed horizontal course, together with their marked time- 
duration ratios that require more attention, as revealed from the numerical studies by Liu, Ellington et 
al (1998), and by Sun k, Tang (2002) on 3-D N-S equations. This last issue on stroke asymmetry has 
been further pursued by Yu & Tong (2005), attaining firm evidences that the aerodynamic asymmetry 
between the down- and up-strokes is a key to understanding the fluid physics of generation of the lift 
and thrust in insect flight. In addition, Yu, Tong & Ma (2005) showed how the unsteady high lift can 
be enhanced by raising the added inertial effects, suppressing the leading-edge flow separation, and 
accelerating the trailing-edge vortex separation. Further, Bao et al.(2006) have studied the effect due to 
passive wing flexibility, discovering it being another mechanism for improving insect flight performance. 
On this issue, Yang et al.(2008) have pursued the idea of integrating unsteady fluid mechanics and 
deforming body dynamics as a cross-discipline between fluids and viscoelastic solids. 

4.3. Computational Fluid Dynamics (CFD) of insect flight 

The problem of insect flight can be resolved numerically by CFD. An integrated model with multi- 
blocked panel grid geometry has been constructed by Liu (2009) to evaluate inertial and aerodynamic 
forces, torques and powers expended for flapping flight in the meso Reynolds number range, with case 
studies demonstrating its accuracy. The effects due to corrugation of insect wings have been found by 
Luo & Sun (2005) to be small on the aerodynamic lift of fruit fly, honeybee, dragonfly, etc. ten kinds 
of insects, and that the effects of aspect ratio (here 2.8 < AR < 5.8) are also surprisingly small. A 
simple model of LEV is given by Maxworthy (2007) to explain that a spanwisc pressure gradient and 
variations in the circulation are responsible to forming a downward propagating vortex ring. The 3D 
flow past a low AR wing translating impulsively for 300 < Re < 500 has been studied numerically by 
Taira & Colonius (2009) for the effects of aspect-ratio, incidence angle, and planeform shape on the 
wake vortices and force generation, with the results validated by oil-tank experiments. 

Related to hovering of insects, PuUin &; Wang (2004) investigated the aerodynamic forces on a 
flat plate accelerating from rest at flxed incidence in a 2D power-law flow using a new analytically and 
numerically related method so specially developed. Their results brought forth a mechanism underlying 
dynamic stall based on a combination of unsteady vortex lift and added mass, under a specific set of 
premises. The force per unit span predicted by the vortex theory is evaluated for parameters typical of 
insect wings and is found to be in agreement with numerical simulations. Estimates for the shed vortices 
and the size of the starting vortices are also obtained, with the significance of the results deliberated as a 
mechanism for insect hovering flight. More recently, the hovering of insects has been simulated by Bos et 
al. (2008) who applied a 2D Navier-Stokes solver to fruit fly hovering at Re = 110 for comparisons with 
other models on the lift and drag, with results revealing that the forces based on their wing kinematics 
differ notably from those based on other simplifled wing kinematics models. 

4.4. Experimental Investigations on insect flight 

Experimental flow visualizations and aerodynamic force calculations of three types of clap-fling 
mechanism have been made by Sohn & Chang (2007). Their tests cover the 'fiing-clap-pause' and 'clap- 
fling-pause' two types, and their corresponding computations applied a 2D Navier-Stokes equations 
solver, flnding good agreement between the two results for the main features about the two families of 
vortex pairs generated and moving off. Results show that the leading-edge vortex bubbles are dominant 
features in causing a large negative pressure region near the leading edge, hence the high lift effects. For 
two-pair wings of large aspect-ratio such as the dragonfly, their flight was examined by Levy and Seifert 
(2009) for i?e < 8 X 10^, finding that the LEV produces relatively weak wake vorticity, thereby reducing 
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the drag and increasing the flight performance, as is well attested by wind and water tunnel tests. On 
insect hovering, Lu, Shcn and Su (2007) made an electromechanical model for wind-tunnel tests, finding 
that (1) spanwise flow is conspicuous, (2) the main leading-edge vortex is led by a secondary dual one 
of the same sense, (3) development of the LEV is somewhat hampered by wing translation, and (4) the 
interaction between fore and hind wing pairs is detrimental to LEV development. 

5. FISH SWIMMING AND BIRD FLIGHT IN NATURE - OPTIMUM CONTROL 

So far, the aquatic and aerial animal locomotion has been addressed with the fluid media primarily 
at rest. In wild nature with the prevailing fluids richly endowed with intrinsic energy, such as in tumbling 
water waves, torrential rivers, or in whirling winds, it should be of interest to explore how these animals 
take advantage of such energy resources by applying their talents evolved over agelong times. Intuitive 
observations often suggest that they are deft players in this game. The famed view of an albatross riding 
on sheared wind {U{z), 0,0), with dU/dz > 0, gaining for its body a velocity {u',v',w') in the frame of 
reference moving with the wind, but without doing any work. This is accomplished mainly by taking 
the advantage of the work done by the inertia force at the rate of —mw'dU/dz, m being its mass, as 
explained in essence by Lighthill (1974). Also, the albatross and pelicans are seen to skim ocean waves 
over long distance with only gentle twisting of their wings, as expounded for the underlying mechanisms 
in detail by D.M. Wu et al.(2005). Porpoises and dolphins have also been watched to ride on bow waves 
in front of a cruising ship with only their tail fluke aptly bent. 

5.1. Migratory Pacific salmon 

Another intriguing story of the Pacific Salmon returning from life-long living in the Ocean to 
swim thousand miles up the Columbia River to the shallow brooklets of their birth place to spawn was 
investigated by Osborne (1960). He tagged some migrating salmon at the Columbia Delta, and later 
recaptured them from fish ladders at hundred miles intervals up the river to acquire data on the time 
of travel and loss of body fat. These salmon are known to stop feeding once having entered fresh water, 
relying on the stored fat as the sole energy supply. He found, quite surprisingly, that a swollen river 
due to a torrential rain upstream did not slow those salmon down by that much of a margin (for known 
metabolic energy they just expended) according to frictional resistance of water being proportional to 
the square of the stream velocity relative to fish. Several possible explanations were suggested, including 
the prospects that the salmon would be able to extract the energy of turbulent eddies for conversion 
into thrTist, navigating through slower stream route, and other ways and means. 

5.2. Extraction of flow energy 

Such fascinating discoveries have inspired further studies on extraction of fiow energy in general. 
Earlier, the interest was first related to resolving problems of fiutter of airplane wings fiying in gusty 
air. This requires critical consideration on the conditions of extracting energy from the stream, either 
actively or passively. In the case of flutter in uniform stream, it is usually assumed that the engine 
maintains the flight speed constant regardless of the flutter-created inertial drag. In waving streams, 
however, the flutter may become self-excited as studied by Kussner (1935) and Garrick (1936). 

5.3. Emulating applications 

Generalizing wings flying in air to hydrofoils cruising in water, this important extension was led by 
Georg Weinblum (1954), a mentoring master in naval hydrodynamics, who developed an approximate 
theory for a hydrofoil heaving and pitching in regular water waves of small amplitude. His views and 
conceptions, broad and profound as usual, offered foresights set on the instability and control of the 
naval vehicle, to both of which a clear understanding of the energy balance is of great importance. 
Further studies on various related issues include that on optimum swimming shape functions by Wu 
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(1971a,b,d), on extracting intrinsic flow energy by Wu (1972), on animals in wavy streams by Wu &: 
Chwang (1975), on utilizing unsteady winds and ocean currents by Wu (1980) and others. In brief, the 
central idea is to find the optimum wing motion with respect to the intrinsic flow field by maintaining the 
feathering parameter nearly at its optimum on control theory. For instance, there have been literatures 
on a family of migratory geese flying in a wedge form, led by a strong member, followed by others on 
the outer upwash sides. Similarly, a school of fish also swim in diamond patterns by taking advantage 
of the vortices with forcwash trails as so suggested in Fig. 2 and in Lighthill (1973a, Fig. 12). 

6. ENERGETICS OF FISH LOCOMOTION 

Animal locomotion is a multi-discipline in which the frontiers of mechanics, biology and related 
fields converge. It benefits from scientists working in intimate collaboration, with keen interests in 
continuing communication and learning from one another. Such collaborative studies can be surprisingly 
rewarding, as exemplified here especially on energetics of animal locomotion. 

After the first International Symposium on 'Swimming and Flying in Nature' held at Caltech in 
1974, another international symposium was organized by Lighthill and Weis-Fogh, held in September 
1975 at Cambridge University on 'Scale Effects in Animal Locomotion,' a new frontier of great impor- 
tance. There, by invitation, Wu (1976) took up a new study on the scaling of aquatic animal locomotion, 
while emulating the comprehensive paper by Lighthill (1976) on the scaling of aerial locomotion. It is 
of interest to give a recount of the proceedings then and its sequel. 

For steady fish locomotion, a mean thrust T is required of the fish to just balance the viscous drag 
D of water on fish, by Newton's law, in sustaining a swimming velocity V, 

T = D (force balance). (49) 

Here, fluid mechanics is suitable in attaining thrust T, by analysis and computation on reactive theory 
to desired accuracy as expounded in §2. For viscous drag D, physiological schemes are suitable for 
evaluating the metabolic conversion from oxygen consumption by fish to muscular power with accurate 
results. So, close collaborations on this joint task between mechanics and biology should be evidently 
fruitful, as is actually a case of splendid success. Otherwise, it is well known that various attempts 
in measuring viscous drag D of immobile fish or their artifical models in wind or water tunnels have 
produced data so much scattered as being useless (for reviews see Gray 1968; Webb 1975). 

To proceed on the joint approach, it begins with the principle of energy conservation asserting 
DV = P-E 

— TjmPm (conservation of energy) , (50) 

signifying that the animal's mean rate of working against drag D is provided originally by the power 
P expended by the fish in performing the work at swimming speed V while casting off, irreversibly, 
some wasted energy E to the fluid, with this work related to the measured net biochemical metabolic 
power P^ consumed all for swimming (after subtracting off the basal metabolic rate which is unrelated 
to swimming movement), in terms of rjm for the 'muscle efficiency' by which the biochemical energy is 
converted to muscular power for active swimming. 

6.1. Metabolic rate and scale effects 

Here, it is important to point out that the metabolic rate in fish locomotion depends on several 
key factors including (a) level of activity, (b) temperature of the water (affecting the concentration of 
dissolved oxygen), and (c) such factors as pre-conditioning (a period of fasting and exercise prior to 
the test to separate energy spent in digesting and interior physiological processes), as well as the state 
of maturity and gender of the specimen. The level of activity is generally hard to characterize. For 
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the Pacific Sahnon, Brett (1963, 1965a, b) has found three distinct levels of performance: the sustained 
(speeds that can be maintained almost indefinitely), the prolonged (speeds maintained for 1 to 2 hours 
with a steady effort, but leading to fatigue), and the burst (speeds achieved at maximum effort lasting 
for only 30 seconds). The physiological basis of each of these levels is different and further depends on 
the scale of body size, such that the measured metabolic power assumes the general formula: 

P^ = am\ b = 0.775, 0.846, 0.890, 0.926, 0.970, (51) 

where m is body mass, proportional to the cube of body length i, and a, b are constant indices, with 
b reported by Brett for the salmon Oncorhynchus nerka, all with an error bar of ±0.145 or less for 
the standard, 1/4- max, 1/2-max, 3/4- max and Vmax-speed, respectively. Here Vmax denotes the speed 
sustainable for 60 min in fresh water at 15°C and the standard state is obtained by extrapolating the 
O2- consumption versus velocity curve to zero swimming speed. Here we note that these b values are 
all greater than the "surface law" (6 = 2/3) and increases with activity level to approach direct mass 
proportionality (6 = 1). Brett's endeavor has made contributions of lasting value to render studies on 
scaling of fish locomotion more systematic and comprehensive, as we now proceed to show. 

6.2. Scaling of svi^imming velocity and energy cost 

Expressing the drag D in terms of the drag coefficient Cn, 

D = ^pV^SCd, Cd = CniRe), Re = Vl/u (52) 

for fish with surface area S (proportional to its length i squared) in water of density p and kinematic 
viscosity {u = 0.0114 cm-^s~^ at 15°C), with Cd dependent on the Reynolds number Re and body 
configuration. From the energy balance (|50p . metabolic scaling (|5ip . and (|52p . it follows that 

For the carangiform and lunate-tail modes of swimming, both modes being well-streamlined, Wu (1976) 
assumes that Cd may be approximated by that for steady flow over a flat plate at equal Reynolds 
number, so that Cd is proportional to i?e~^/^ or Re~^^^ according as the boundary layer is laminar or 
turbulent, respectively, hence yielding the scaling law: 

y = const.^'^, with /3 = ^(26-l), ^(56-3), 6-^, (54) 

according as the boundary layer is laminar, turbulent, or with Cd = const., respectively. The last case 
(of quite large but constant C^i) is for the hypothetical case when flow separation occurs in the cross 
flow past an undulating fish body in a braking-type maneuver. 

In addition to providing the detailed metabolic rate of a single size group, a comparative phys- 
iomechanical study has been made possible by Brett's data (1965) on P and swimming speed V for 
sockeye salmon of five different size groups. The scale effects are exhibited in the Cp ^ Re relationship 
as shown in Figure 8 with five size groups (in dashed lines), each covering the four activity levels as 
marked. The least-square-error fit for each activity level of the five size groups brings out the scale 
effects. In this log-log plot of the Cp ~ Re lines, the slope reads -0.61, -0.54, -0.49, -0.41 for 1/4-max, 
1/2-max, 3/4-max, and full level of activity, respectively, for fish over five length scales. The result of 
this analysis is already an important discovery because the 1/2-max and 3/4-max level lines are situated 
across that of the reference Cd ~ Re line for the laminar (with slope —0.50) in sharp contrast with the 
turbulent Cd (with slope —0.20) of a smooth flat plate shown in lower Figure 8. Also shown are the 
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Figure 8: Scale efTects in the measured metabolic power coefficient Cp for swimming exhibit variations with Reynolds 
numbers Re for five different size groups of sockeye salmon: (solid lines): least-square-error fit to Cp at specific activity 
levels; (dashed lines): Cp variations with various activity levels for each size group; (dash-dot lines): theoretical thrust 
coefficient Ct - Figure adapted with Brett (f965a,b); courtesy of Academic Press. 



thrust coefficient Ct lines derived with using the observed data of Pyatetskiy (1970) and Webb (1975). 
The slopes of these Ct lines are seen to be qualitatively similar to those of the estimated Cp curves for 
fish of each scale group (see Fig. 8 ). This scrutinizing scaling therefore brings forth the conspicuous 
result that the viscous drag on all these specimens is primarily associated with laminar boundary-layer 
flows as so attested by this solid and sound evidence. 

In addition, we introduce the ^specific energy cosf, £, defined for fish of weight mg as 

a dimensionless parameter which provides a valuable measure of the relative merit of the propulsive 
system, signifying the energy expenditure in transporting unit weight over unit distance. It was in- 
troduced by Karman and Gabrielli (1950) in evaluating the comparative merits of fourteen classes of 
transportation vehicles and animal locomotion, and used by Schmidt-Nielsen (1972), Tucker (1975) and 
others in studies of comparative physiology, and also by Weihs (1973) in his study. 

These results for the scaling of metabolic rate and swimming velocity may now be combined to 
evaluate the scaling of the specific energy cost. From (jSip . (j54p and (|55p we find that 

£^ = P^/mgV = const. m-'^, 7=1-6-^/3/3; 7 = 0.32,0.29,0.25,0.21, (56) 

where the values for 7 are deduced by using the experimental data in Fig. 8 and the values of b in 
(51) for the activity level at standard, 1/4-max, 1/2-max, and 3/4-max level of activity, respectively, 
showing how the specific energy cost £ varies with body mass accordingly. These values of 7 and the 
observed mean of /3 = 0.5 are plotted versus b (as a measure of the level of activity) in Figure 9. In 
these experimental data, the basal metabolic rates have been subtracted off (cf. Wu & Yates 1978). 
Also shown are the similarity predictions of (3 and 7 evaluated for the three distinct reference states 
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Figure 9: Comparison of the observed values of the scale similitude /3-line and the specific energy cost 7-line (adapted 
with Brett 1965a, b) with the predictions based on three reference states for laminar or turbulent boundary layer, and 
Cd =constant for the braking case of separated cross flows. The b values have the basal metabolic rates subtracted off, 
and the bio-chemically observed /3-line and 7-line both intersect only with their corresponding laminar boundary-layer 
reference lines, and both precisely at the 1/2-max activity level, with the intersection points circled in red. 



of laminar, turbulent and separated flows. Here, a comparison between the experimental data and the 
similitude calculations brings forth far-reaching results. First, the bio-chemically observed f3-line and 
"f-line both intersect only with their corresponding laminar reference lines of the three flow reference 
states, and both at exactly the 1/2-max activity level, with the intersection points circled in red. This 
therefore shows that the boundary layers of these swimming fishes should have to be evidently laminar at 
these high Reynolds numbers (8 x 10^ < Re < 1(f), and further shows that the specimens tested would 
have to be swimming about 1/2-max of their maximum sustainable level. 

7. CONCLUSION 

The present expository review has delineated classical studies in five major areas of animal loco- 
motion: fish swimming, bird flight, insect flight, propulsive motion control, and locomotion energetics. 
Each of these areas has marked with its own nature, outstanding properties, and significant advances. 
For fish swimming, passing an undulatory wave of speed c distally down its flexible slender body to 
swim at velocity U slightly less than c is the primary key feature. To the top performers, the appended 
dorsal, ventral and caudal fins evolved in intricate ways display challenging nonlinear and interactive 
mechanisms that can now be expounded by theory and experiments developed. Similar statements can 
be properly applied to bird flight. The attractive flapping of bird wings stretched straight in downward 
pronation stroke and the bending and twisted wings in uplifting supination stroke can now be well 
modeled mathematically for analysis. Fully nonlinear and unsteady theory can be applied for evalua- 
tion of the aerodynamics of a flexible wing moving in arbitrary manner along arbitrary trajectory for 
further applications. The subject on insect flight has appeared impressive in the diversity of countless 
species, each with differing mode of flapping wings of unique shape and beating frequency. The relevant 
aerodynamic category falls in the range of meso Reynolds numbers of 10^ < Re < 10^. However, the 
viscous effects are rendered negligible by the high values of the Strouhal number resulting from the 
high frequency of the insect wing flapping. In addition, recent studies in this field have discovered 
an outstanding characteristic mechanism rather universal in the insect world. And this marked the 
wing flapping at finite incidence angles of 20° — 40°, causing flow separation from both the leading and 



27 



trailing edges in forming a conspicuous leading-edge vortex bubble to generate unsteady high lift. This 
mechanism has been adopted to make mathematical models for insect flight studies. 

On the energetics of animal locomotion, the original idea is brining the efficient and effective part of 
mechanics (in evaluating thrust delivered by fish) and of biology (capable of metabolic measure of mus- 
cular power for counter-balancing drag on swimming fish) for unified studies. The results are rewarded 
by the invaluable evidences achieved, which would be utterly impossible otherwise. In concluding §6, 
it has been shown, with Fig. 9, that a comparison between the experimental data and the similitude 
predictions brings forth far-reaching results. First, the bio- chemically observed f5-line and ^-line both 
intersect only with their respectively corresponding laminar reference lines of the three flow reference 
states, and both intersection points fall precisely at the 1/2-max activity level. These results therefore 
furnish the valid and sound evidences showing resoundingly that the viscous drag on fish swimming in 
nature is associated primarily with laminar boundary-layer flows. And they are thus providing a valid 
and rigorous resolution to the 'fish swim paradox' proclaimed earlier by Sir James Gray. That para- 
dox has been a mechanics-biology joint stimulus that has attracted much of the diligent and dedicated 
studies over the decades, now bearing a fruit to both mechanics and biology disciplines for celebration. 
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